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TRANSFER FUNCTION OF A 0T16D NOZZLE 
DUE TO ACOUSTIC AND CONVECTED ENTROPY FLUCTUATIONS 

by 

J. H. Miles 

National Aeronautics ana Space Administration 
Lewis Research Center 
Cleveland, Ohio 
ABSTRACT 

An acoustic transmission matrix analysis of sound propagation in a vari- 
able area duct with and without flow is extenoed to include convected entropy 
fluctuations. The boundary conditions used in the analysis are a transfer 
function relating entropy and pressure at the nozzle inlet and the nozzle exit 
impedance. The nozzle pressure transfer function calculated with this analy- 
sis is compared with JTlbU turbofan engine nozzle data. The one dimensional 
theory for sound propagation in a variable area nozzle with flow but without 
convected entropy is good at the low engine speeds where the nozzle exit Mach 
number is low (M=0.2) and the duct exit impedance model is good. The effect of 
convected entropy appears to be so negligible that it is obscured by thy 
inaccuracy of the nozzle exit impedance model, the lack of information on the 
magnitude of the convected entropy and its phase relationship with the pres- 
sure, and the scatter in the data. It appears an improved duct exit impedance 
model is required at the higher engine speeds where the nozzle exit Mach num- 
ber is high (M=0,56) and at low frequencies (below 120 Hz). Nevertheless, the 
results show the one dimensional sound propagation model does apply to the 
0T15D turbofan engine nozzle data. 


NOMENCLATAURE 


matrix 

0 

A area, m 

Cq isentropic speed of sound, m/s 

c„ gas specific heat at constant pressure, O/kg K 

r 

Cy gas specific heat at constant volume, J/kg K 

f frequency, Hz 

i (- 1 )^^^ 

k wave number, u/c, m"^ 

L duct length, m 

Mq Mach number 

(MW) gas molecular weight 

111 d ln(A)/dx, iif^ 

P transformation matrix 

p pressure, Pa 

R gas constant, J/kg K 

9 

S temperature power spectral density, K /Hz 

s entropy, J/kg K 

T transmission matrix 

t time, s 

u velocity, m/s 

x Cartesian coordinate, m 

Y acoustic state vector 

1 acoustic impedance, mks Rayles 

Y gas specific heat ratio 

AX region length, m 

X eigenvalue 
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p gas density, kg/ni 

0 temperature, K 

u) angular frequency, radians/s 

Superscripts 'and subscripts; 

0 referg“nce state quantity 

1 perturbed quantity 

INTRODUCTION 

Internally generated turbojet engine noise is important at low power 
settings where the core noise level usually excqeas that of the pure jet 
noise. This noise can be attributed to several sources, such as combustion, 
turbulence in the turbine or nozzle, and flow incidence upon struts. 

The combustion process produces two noise generation effects. The first 
can be attributed to the fluctuating heat release rate which produces isen- 
tropic pressure waves that propagate through the engine and are radiated to 
the far field. The second can be attributed to local regions of hot gas 
(entropy variations) which convect through the engine with the mean flow. The 
entropy variations are not a direct noise source since they have no associated 
pressure disturbance. However, as the entropy convects through regions with 
area variations acoustic waves are generated. Furthermore, the isentropic 
pressure waves generated by the combustion process also interact with area 
variations. The resulting waves represent a combination of waves due to both 
sources of change. 

As part of a core noise research program conducted at the NASA Lewis 
Research Center, far field and internal pressure measurements were made near 
the entrance and exit of a JT15D core exit nozzle over a range of engine power 
settings (Ref. 1). Using these measurements nozzle pressure transfer 
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tunctionj. were calculated. Similar aata tor the YF102 turbotan engine were 
available. However, the core nozzle exit and inlet areas of the YF102 are 
nearly equal while the exit area of the JUSD core nozzle is half the inlet 
area. Since an entropy perturbation must convect through an area channge to 
produce an acoustic pressure, only the 0T15D data are used herein. 

In this report, first the acoustic transmission matrix approach developed 
in Ref. 2 to study isentropic sound propagation through a variable area duct 
or nozzle carrying a compressible subsonic flow is extended to induce con- 
vected entropy effects. Then this analytical propagation model is used to 
calculate a nozzle exit pressure transfer function. Last, the measured and 
calculated nozzle pressure transfer functions are compared. 

BACKGROUND 

Entropy noise is one of the many sources of excess noise reviewed by 
Ffowcs Williams in Ref. 3. Many studies of noise due to entropy convecting 
through area changes have been made. The problem of noise production due to 
convected entropy fluctuations in a choked nozzle was studied theoretically by 
Candel (Ref. 4) and Marble and Candel (Ref. 5). Noise production by this 
mechanism in a subsonic nozzle was studied theoretically by Bohn (Ref. 6), Lu 
(Ref. 7), and Ffowcs Williams and Howe (Ref. 8). Entropy noise has been 
studied experimentally by Zukoski and Auerback (Ref. 9), Auerback (Ref. 10), 
Bohn (Ref. 11), Strahle, et al. (Ref. 12), and Muthukrishnan, et al. (Ref. 13). 
A brief but comprehensive review of the early work is given in Kef. 5. 

Noise production by entropy fluctuations passing through turbines of cur- 
rent aircraft engines has been studied by Pickett (Ref. 14) and Cumpsty and 
Marble (Refs. 15 and 16). 

The governing equations applied herein were originally used by Tsien (Ref. 
17) and Crocco and Cheng (Ref. 18) to study combustion instability in a rocket 
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combustion chamber. The approach taken herein Is to decompose the entropy 
inhomogeneity into a Fourier distribution of entropy variations. These 
equations and this approach also were used in Refs. 4 tc 7. 

In the study of noise due to convected entropy fluctuations, Candel (Ref. 
4), Candel and Marble (Ref. 5) and Bohm (Ref. 6) introduce the assumption of 
linear mean velocity distributions in the nozzle used by Tsien (Ref. 17) which 
allows some analytical simplifications. The theoretical study of this problem 
was simplified by Lu (Ref. 7) when he used an exponential nozzle in formulat- 
ing his entropy noise model. The analysis presented herein is more general 
since it is applicable to arbitrary nozzle shapes and velocity profiles. 
Furthermore, it is simple to use and yields results directly comparable with 
measurements. 

In this paper, the governing equations are derived by linearizaton and put 
into state variable form using acoustic pressure, particle velocity, and en- 
tropy as the state variables. Then, in order to solve the resulting differen- 
tial equations, the nozzle is divided into regions in which the coefficients 
of the differential equation can be made constant by assuming an exponential 
area variation. A transmission matrix is obtained for each region by solving 
the acoustic state variable differential equation using the eigenvector matrix 
exponential method described by Ogata (Ref, 19). Other procedures for comput- 
ing the exponential of a matrix are reviewed in Ref. 20. The transmission 
matrix of the nozzle is the product of the transmission martices for each 
segment. This procedure is used in Ref. 2 for the case of zero convected 
entropy. It is experimentally verified for the case of sound propagation in a 
variable area duct without flow in Ref. 21. 

In order to calculate the nozzle pressure transfer function, two boundary 
conditions are needed. For one boundary condition, it is assumed that the 
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nozzle exit impedance is given by a model proposed by Lumsdaine and Ragab 
{Ret. 22). The second boundary condition is the ratio ot the convected 
entropy perturbation to the acoustic pressure at the nozzle exit as a function 
of frequency and it is determined empirically. 

THEORY 

The following assumptions are made to obtain the plane wave acoustiv: 
equations: the duct has an area profile A(x); the acoustics are described by 

first order, linearized equations of continuity, momentum, energy and state 
for isentropic sound propagation. 

State variable equation 

The one dimensional flow in the nozzle of a nonviscous, nonconducting, 
perfect gas is oescribed by the familiar equations of continuity, momentum, 
energy and state (Eqs. (1) to (4) of Ref. 2). 

The fluctuations in the duct are assumed so small that the flow in the 
duct is only slighthly perturbed. Consequently, the instantaneous quantities 
can be written in terms of a unperturbed stationary flow quantity designated 
by 0 and a small perturbed quantity designated by 1 (Eqs. (5) to (8) in 
Ref. 2). 

Substituting the equations for the instantaneous quantities into the 
equations of continuity, momentum, and energy yields at zeroth order one 
system of equations (Eqs. (9) to (11) of Ref. 2) and at first order another 
system of equations (Eqs. (12) to (14) of Ref. 2). 

For the system under consideration the gas is assumed to be a perfect gas. 
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The undisturbed flow is assumed to be steady and isentropic so that s^ 
and (p/pp are constant in time. Consequently, the energy perturbation 
equation (Eq. (14) of Ref. 2 ) can be integrated to 


0 (t,x) n {t,x)/Cy 


Pi (t,x) 

“p^W* 


Y Pi (t>x) 

“••TTixr” 



dx 

e^lxT. 


( 1 ) 


Thus an entropy fluctuation is convected with the gas flow and its value 
remains unchanged from a value o„(t) determined at the nozzle entrance, x 
= 0.0 . 

In Ref. 2 the entrance entropy fluctuation was assumed to be negligible. 
Consequently, the density perturbaton was related to a pressure perturbation by 


Pi - Pi/Cq (2) 

However, in this paper, the entrance entropy fluctuation is assumed 
significant and 


Pi Po h 

•"I T ” T 

^0 S % 


(3) 
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Thfi perturbed state variables are the acoustic pressure, the particle 
velocity, and the entropy. These quantities are assumed to vary with time as 
exp(-iu)t). Substituting Eq. (3) Into the Fourier transform of Eqs. (13) to 
(14) of Ref. 2 yields the following equation in matrix form: 



The differential equations used in Refs. 4 and 5 to model a variable area 
duct with acoustic and convected entropy perturbations differ from the 
familiar first-order acoustic equations in two respects. They have additional 
terms on the left hand side which represent the changes in the field due to 
gradients in the duct geometry, mean flow and density. Furthermore, they have 
dipole source terms on the right hand side which arise from the interaction of 
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the entropy fluctuation with gradients of the duct geometry and mean flow. 
Consequently, in these analytical models the entropy fluctuation is viewed as 
a source term and the equations are treated as being nonhomogenous. Thus, in 
order to solve these equations an acoustic impedance must be specified at the 
start and end of the source regie", and the entropy fluctuation must be 
specified at a point in the region. 

However, in this paper the first-order equations are written with the 
terms which arise from the interaction of the entropy fluctuation with the 
gradients of the duct geometry and mean flow on the left hand side. As a 
result the equations have no source terms. Thus only propagation occurs and 
the equations are homogeneous. Thus, in order to solve these equations only 
two boundary conditions at an end of the propagaton region are necessary. As 
a result, specifying an acoustic impedance boundary condition and the ratio of 
the convected entropy perturbation to the acoustic pressure as a function of 
frequency at the same end of the propagation region is sufficient to generate 
a solution. 

The gradients of the logarithmical variation of velocity and density are 
related to the logarithmical variation of area as follows; 


d in(u^) 
dx 




d 1n(A) 
dx 


d In(pQ) 

3x 



1n(A) 

dx 


( 5 ) 

( 6 ) 
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Substituting Eqs. (b) and (6) into Eq. (4) and solving for dY/dx where Y n 
(Pl> s^) yields 




dY/dx e BY 

(7) 

where 

the elements of the B matrix are: 


^ 11 “ 

Mo (-iko) 

4 . .. '^P d ln(A). 

( 8 ) 


1 



j. '’o'^o ^ Mo 0 in(A) 

(9) 

1 - Mo 

,, 2 dx 

^13 = 

p 

^0 0 

'^o d ln(A) 

( 10 ) 

*^p 1 

- “ 
0 

B^l 



( 11 ) 

^ 0 % - 

M^) 

&22 = 

Mo (-iko) 

d ln(A) 

( 12 ) 

1 - 

0 

« 2 dx 

(l-M^ 

^23 = 

0 


(13) 

B 3 I = 

0 


(14) 

B 32 = 

0 


(15) 

B 33 " 

i 0 )/U^ 


(16) 
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where d w/e^. Note that the 3x3 B matrix in Eq, 7 differs from the 
^x?. acoustic propagation B matrix discussed in Ref. 2 by the addition of 
matrix 'element ^13 ViJ <ch gives the change in the gradient of the acoustic 
pressure due to the convected entropy and matrix element B^g which gives 
the change in the entropy gradient due to the convected entropy. The acoustic- 
convected entropy state variable differential equation is given by Eq. (7). 

This equation will be solved next. 

Method of solution 

In the previous section a general acoustie-convected entropy state vari- 
able differential equation was derived (Eq. (7)), In this section an approx- 
imate solution is found, A method of solution similar to the one used in this 
study is presented in Ref. 2, However, in Ref. 2 it is applied only to the 
acoustic equations. For clarity, it is briefly discussed and those modifi- 
cations nedessary to adopt it to these acoustic-convected entropy equations 
are given. The first step in obtaining a solution to Eq. (7) is to divide the 
duct into a number of regions or subsections. The region size is selected so 
that the area variation in each region can be approximated by an exponential 
area variation. Consequently, in a given region 

= constant 
and 

A = e"“ 


(17) 

(18) 
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With this assumption the va'<a of d ln(A)/dx is a constant in the region 
given by 

m o a ln(A)/dx (19) 

The B matrix in Eq, (7) is now independent of x in the if'egion where it is 
to be evaluated. Thus, the solution to Eq, (7) for the jth region is 


Y (Xj+ ax) = (e^ ^^)^ Y (x^) 


( 20 ) 


where x is the length of the jth region and exp(B ax) is known as the matrix 
exponential. The matrix exp(B ax) is the transmission matrix of the jth 
region, (T)j, so that 



S .£=0 


11 

'•^12 

Ti 3 

21 

”^22 

”^23 

i 

0 

^33 


( 21 ) 


Consequently, the transmission matrix T for the variable area duct is found 
from 

N 

Y (x=L) - TT 
j = 1 


where 


N 

I - S (23) 

j = 1 ^ 

and N is the number of nozzle regions. 
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Evaluation of matrix exponential 

While the series given in Eq. (21) defines the matrix exponential, it was 
not used for computer calculations since the convergence of the series may be 
slow. The matrix exponential was evaluated from 


where 


= C“^ Pe^ P”^ C 


D = P“^ ^P 
y= CBC“^ 


1 


C = 



0 


0 


= 

0 

0 





(24) 


(2b) 

(26) 


(27) 


and P is determined below such that the matrix D is diagonal. The P 
matrix then consists of eigenvectors of . 

The matrices in Eq. (24) as used in Ref. 2 are 2x2 since acoustic pressure 
and particle velocity are the only variables. However, as used herein the 
matrices are 3X3 since entropy perturbations are being included. Consequently, 
the C matrix contains a new element. Furthermore, the eigenfunction P 
matrix used herein has a new definition dependent on the eigenvalues 
Xp X 2 > and Xg of the matrix. 
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For j G 1 and 2 the column eigenvector of P is selected as follows; 





(28) 


^32 “ *^21 


(29) 


P. 3 GO 


(30) 

However, the third column of P is selected to be 




P 3 I = “ y ^3 


(31) 


P 32 = ^^21 ^<^13 


(32) 


^33 " (^^22"^-3^ “ '^IZ 

s^Z\ 

(33) 

The values of 

Xp X 2 » ^3 eigenvalues of 

the matrix. 


The values of 

and x^ ane given by the roots of 




2 

X + bx + c = 0 


(34) 

where 





b 


(36) 

and 





c = s4ii s^22~ ^IZ ^Zl 


(36) 

The value of X 3 

is given by 




^3 = '“/‘'o 


(37) 
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It can be verified that 


hence 


0 Ci 



c 




(38) 


(39) 


Nozzle exit acoustic boundary conditions 
The nozzle calculations shown herein are made using the following duct 
exit impedance 


p c 
^0 0 


= Mq + 1 - 2 (v)/v - i (v)/v 


(40) 


where 


V 


2 

0 d 



(41) 


(v) is the Bessel function of the first order and first kind, (v) 
is the Struve function of the first kind and first order, and r^ is the 
nozzle exit radius. The duct exit impedance was derived in Ref. 22 for a 
circular duct with flow having an open end fitted with an infinite acoustic- 
ally rigid flang. Note that in the case of outflow, M is replaced by 
- M in Eq. (40). 
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Calculation of nozzle transfer function 
Tlie transfer function between point (2) near the core nozzle entrance and 
point (1) near the core nozzle exit is calculated using the acoustic pressure 
at the core nozzle entrance, p^ (x=0.0), the acoustic impedance at the core 
nozzle entrance, Z(x=0.0), and the convected entropy at the core nozzle, 

Zg (x=0.0), as follows: 


P (2) 


and 


P (1) 


■[ 

■[ 


. , , ^12 ( 2 ) 
hi 1 rx=o) 


h? (2) 

hi ■* I (xiDTO) 


Pj (x=0.0) + lj 3 (2) Ej{x=0.0) 


Pj (x=0.0) + T33 (1) 2j(x= 0.0) 


(42) 


(43) 


The transfer function, H(f), is the ratio of p(l) to p(2) 


H (f) = p (l)/p (2) 


(44) 


Consequently, the major effect of convected entropy depends on the transfer 
function relating entropy and pressure at the nozzle inlet. 

However, the acoustic impedance at the core nozzle entrance is also a 
function of the entropy. The acoustic impedance at the exit is assumed to be 
independent of the convected entropy. The acoustic impedance at the core 
nozzle entrance Z(x=l) is found from the acoustic impedance at the core 
nozzle exit, the convected entropy at the inlet and the convection time oelay 
using the following relationship: 
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The absolute value of (xcO.O) does not determine the nozzle transfer 
function. The value of the convected entropy and the convected entropy-pressure 
transfer function at the entrance determine the nozzle transfer function. The 
value of Pg (xcO.O) is selected so that for zero convected entropy the 
calculated value of the auto pressure spectra at points one and two are similar 
the the measured pressure auto-spectra at these points. Since the measured 
pressure auto-spectra may be due to both acoustic and nonacoustic pressure 
perturbations they were not used directly. 

In order to define an acoustic pressure at the core nozzle inlet to use as 

Pg, a fictitious point pressure source was placed at the core nozzle 

1 

entrance. This point pressure source was assumed to have a (f)” spectrum. 

The spectrum level was determined by single constant selected to make the 
calculated and measured auto spectra appear similar. A Green's function 
solution was then constructed in a small region about this hypothetical point 
source using the source region model described in Ref. 23. The upstream 
impedance was constructed by assuming the reflection factor was 0*9. The 
downstream impedance was the core nozzle inlet impedance, Z(x=0.0). 

In these calculations the pressure at the entrance is real. However, the 
phase relationship of the pressure and entropy perturbations is not known. 
Consequently, two arbitrary phase relationships of pressure and entropy were 
used in the computational study. Most calculations were made using 

iu) J dx/u^ (x) 

(x=0) = lEgi e ^ (46) 
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so that the eonvected entropy and pressure were in phase at the core nozzle 
exit. While this relations is unrealistic, it has the advantage that changes 
in the transfer function with changes in the magnitude of could be 
more easily compared to the case where is zero. In order to obtain 
some information on the importance of the phase relationship between pressure 
and flow, some calculations were made using 

£3 (x=0) == lEgi (47) 

For this case the eonvected entropy is real at the entrance and out of phase 

with the pressure at the exit by 03 / dx/u. (x). 

Jq ° 

Various values of the magnitude of the eonvected entropy were used. For 
the cases presented the magnitude of was selected to be large enough 
to show the effects of its presence. 

EXPERIMENTAL DATA 

The Pratt & Whitney 0T15D is a bypass ratio 3.3, two-spool, turbofan 
engine with a rated thrust for 9800 Newtons. The engine core consists of a 
compressor in the form of a 16-bladed impeller, a reverse flow annular com- 
bustor, and a three-stage turbine. The single-stage compressor is driven by 
the high-pressure turbine stage, while the fan is driven directly by the two 
low-pressure turbine stages. The fan has twenty-eight blades and is 
0.534 m in diameter. 
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The test data are discussed in Ref.l. All tests were conducted using an 
outdoor engine test stand. The engine was mounted with the engine centerline 
2.9 m above a hard surface ground plane. The engine was fitted with an inlet 
control device to reduce fan tones by providing a cleaner and more uniform 
inflow than is normally obtained under static test conditions. A photograph 
of the engine mounted on the test stand is shown in Fig. 1. 

Dynamic pressure probes were placed in the engine core at seven different 
locations (Fig. 2). Data used herein were measured using probes at the core 
nozzle entrance and at the core nozzle exit plane. The core nozzle geometry 
and microphone locations are shown in Fig. 3. The transducers used were con- 
ventional 0.635 cm diameter pressure response condenser microphones. To avoid 
direct exposure of the microphone to the severe environment within the core, 
they were mounted outside the engine and the fluctuating pressure in the 
engine core was conmiunicated to the transucers by "semi-infinite’’’ acoustic 
waveguides. These waveguide probes are described in detail in Ref. 24. Data 
were also taken using far-field microphones. Data acquisition and processing 
are described in more detail in Ref. 1. Test conditions for the noise 
measurements used herein are given in Table 1. 

In order to calculate a noise-free acoustic transfer function, the experi- 
mental transfer function is calculated from 

H . . pm .. p* ■ (w 

P(2) p (e) 

where p(2)p (e) is the cross-spectrum between the microphone at the nozzle 

"ic 

inlet and a far-field microphone and p(l)p (e) is the cross-spectrum 
between the turbine exit microphone and a far-field microphone. 
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RESULTS 

Transmission matrix 

Transmission matrices for the duct configuration shown in Fig, 3 with and 
without flow were numerically calculated using Eq. 40. The calculations yield 
an exact solution for the case where the Mach number is zero and an approxi- 
mate solution for the case with flow. The approximation is improved as the 
core nozzle is divided into smaller subsections. The matrices wore calculated 
using 230 equal length subsections. Typical elements of the transmission 
matrices T calculated with and without flow are shown in Fig. 4 where 

Y (x«39.5 cm) = T y (x=0.0) (49) 

For the solution with flow, the zeroth-order system of equations is solved to 
evaluate p, u , p, and e at the start of each subsection using a stagnation 
pressure of 1.0076 atmosphere, a stagnation temperature of 414 C and a mass 
flow rate of 2.73 kg/s. 

The acoustic matrix elements, Tj^^* ^12» ^21’ ^22 

cases with and without flow are similar to the corresponding transmission 
matrix elements of the long exponential nozzle shown in Fig. 2 of Ref. 2. The 
convected entropy matrix elements Tj^2» ^23 ^33 results. 

Matrix element T23 yields no new information for it merely shows that the 
entropy is convected unchanged. However, matrix element Tj^^ shows the 
amount the core nozzle exit pressure is changed by an inlet entropy pertur- 
bation and matrix element T23 shows the amount the core nozzle exit 
acoustic particle velocity is changed by an inlet entropy perturbation if the 
inlet acoustic pressure and particle velocity are zero. Note that the phase 
of the matrix elements Tj^3 and T23 is strongly dependent on the 
convection velocity. 


20 


Transfer function calculations 


The entropy spectrum £_ is related to the gas temperature spectral 

b 

density since 


llil. o e_ . Ikxl 3 i- 
=p «o Po If « 0 


( 50 ) 


where the pressure term is neglected since the acoustic pressure perturbation 
is much smaller than the mean pressure. However, no measurements of the gas 
temperature spectrum in the nozzle were available. The existence of a gas 
temperature spectrum in the nozzle can be inferred from the rms ana power 
spectral density gas temperature measurements in a gas turbine transition duct 
exit presented in Ref. 2b, 

Since no spectral data were available, the acoustic pressure transfer 
function was calculated to show the effect of a particularly large rms 
temperature standard deviation in a single 1 Hz frequency band. The standard 
deviation over a frequency range f^ is related to the temperatue power 
spectral density by 


2 



d f 


(51) 


2 

Consequently, the power spectral density S„ is much lower than e /Hz 

© 

if the standard deviation in the frequency band is 0. 

Most calculations were made using m = 50.0 J/kg K. At the lowest 
test condition when the nominal engine speed is 33 percent of the maximum, 
Cp = 1083 0/kg K, and = 675 K. Consequently, the corresponding rms 
temperature deviation is about 31 K. This corresponds to having an rms 
deviation of 31 K where the temperature is excited over a narrow 1-Hz band- 
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width. If the temperature spectrum was band-limited with noise up to 1200 Hz, 
the rrns deviation of 31 K over this frequency range would imply that the con- 
vected entropy magnitude was 1.44 J/kg K, 

Comparison with experimental data 

Transfer functions were calculated using various models and compared with 
experimental data. Some typical results are presented next. Figure 5 show 
that the transfer function measured at a nominal engine speed of 33 percent of 
the maximum engine speed. Also shown in Fig. 5 are transfer functions 
calculated for the following cases: (1) no flow, (2) with flow but with no 

convected entropy, and (3) with flow and convected ent>"opy with itgi 
= 50.0 J/kg K and the convected entropy and pressure in phase at the nozzle exit. 

The measured transfer function magnitude and phase plots contain much 
scatter. Due to the scatter in the magnitude plot shown in Fig. 5(a), all the 
methods for calculating the transfer function magnitude appear to agree with 
the data over the frequency range from 120 Hz to 720 Hz. Furthermore, all 
three methods fail to yield the correct magnitude above 720 Hz where higher 
mode propagation may be important. Also, only the zero flow case fits the 
magnitude data below 120 Hz. 

The measured transfer function phase angle plots shown in Fig. 5(b) have 
less scatter below 720 Hz. Comparing the calculatea and measured transfer 
function phase angle curves shows that most of the data are bounded between 
the phase angle curve calculated for the flow and no-flow cases for the 
frequency range from 120 Hz to 720 Hz. However, the calculations made for the 
cases with flow are in better agreement. Nevertheless, of the two cases with 
flow, only small differences occur beteeen the cases with and without 
convected entropy. 
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The one-dimensional acoustic propagation theory for a variable area duct 
with flow applied using the duct exit impedance model of Lumsdaine and Ragab 
(Ref. 22) in general produced good agreement with the data. However, the 
results are uot precise enough to indicate the presence or absence of con- 
vected entropy effects. The one-dimensional acoustic propagation theory for a 
variable area duct without flow has been verified experimentally (Ref, 21) and 
it is assumed to be valid with flow. However, the duct exit impedance used is 
only approximate since it was derived for a circular duct with flow having an 
open end fitted with an infinite acoustically-rigid flange. 

Consequently, an improved model for the nozzle exit impedance based on the 
geometry and flow at the OTlbD turbofan engine core nozzle is needed in oroer 
to study the convected entropy effect. 

Figure 6 shows the measured auto-spectra near the inlet and exit of the 
core nozzle at the nominal engine speed of 33 of the maximum engine speed,, 

Also shown in Figure 6 are calculated auto-spectra for the same examples used 
in Fig. 5. Note that the measured and calculated core nozzle inlet 
auto-spectra have a dip near 500 Hz and a peak near 800 Hz, The calculated 
pres'Jure auto-spectra level was adjusted to peak near 135 dB at 850 Hz. 

Figure 7 compares the transfer function measured at 33 percent of the 
maximum engine speed with the transfer function calculated with iz^i » 50.0 
and 100.0 0/kg K and the convected entropy and pressure in phase at the 
nozzle exit. Increasing the level of the convected entropy does not produce 
better agreement. 

The transfer function curves shown in Figs. 5 and 7 were calculated with 
the convected entropy and pressure in phase at the exit. These transfer 
function curves are smooth. In Fig. 8 the transfer function calculated using 
Eq. 47 with the convected entropy and pressure out of phase at the exit are 
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shown* This transfer function is calculated for the case where the conveeted 
rsntropy has a zero phase angle at the entrance. The resulting transfer 
function amplitude and phase angle curves have dips and peaks, Consequently! 
the phase relationship betv)3en the conveeted entropy and pressure must be 
measured in order to know if the dips and peaks in the transfer function phase 
angle data can be attributed to conveeted entropy. 

For the given geometry and the measured auto spectra the value of the 
conveeted entropy had to be 60/kg K in order to see a significant conveeted 
entropy effect. This value might occur in a single frequency 1 Hz band 
because it represents a standard deviation of 31 K in a 1 Hz band. However, 
it cannot exist over a large frequency range or the standard deviation would 
be larger than the average. Consequently, it is more likely that the value of 
the conveeted entropy is of the order 1,44 J/k^^ K over the frequency range 
from 0 to 1200 Hz so that the standard deviation is 31 K in the band from 0 to 
1200 Hz, However, if the conveeted entropy is of this order then no effect 
would be observed. Furthermore, this is true independent of the phase of the 
conveeted entropy. Thus, the effect of conveeted entropy on the nozzle trans- 
fer function at this test condition is so negligible at most frequencies that 
it is obscured by the scatter in the data and use of an incorrect nozzle exit 
impedance. 

These results show that the calculated transfer function with flow and 
with and without conveeted entropy resembles the data measured at 33 percent 
of the maximum engine speed. Moreover, they indicate that conveeted entropy 
could be a factor in defining the transfer function. Furthermore, these 
results indicate that the phase relationship of the conveeted entropy to the 
pressure is as important as magnitude information. 
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Calculations wero also made to determine if the model transfer function 
resembled the data measured at other engine speeds. In Fig. 9 is shown the 
measured transfer function magnitude and phase angle plots at 33 percent, 40 
percent, 55 percent, 87 percent, and 97 percent of the maximum engine speed. 
Also shown in Fig, 9 is the transfer function calculated without and with 
flow, using is^i c 50 J/kg K, and using the convected entropy and 
pressure in phase at the core nozzle exit. Examination of the calculations at 
33 percent of the maximum engine speed shown in Fig. 5 indicates that the the 
transfer function calculated using this value of iZgi is essentially 
identical to that calculated for iZgi equal to zero. Thus, in 
Fig. 9, essentially transfer functions calculated with and without flow are 
compared at the lower engine speeds. 

Examination of the measured transfer function magnitude curves shown in 
Fig. 9 indicates that while the nozzle exit Mach number varies from 0.16 to 
0.56, the magnitude plots do not change significantly over the entire range. 
The measured phase angle curves do show some change with operating conditions. 
The phase angle curves measured for nominal engine speeds of 33 percent, 40 
percent and 55 percent form a similar group. For example, the phase angle 
between 120 and 240 Hz is about 20 degrees for these low engine speeds. The 
phase angle curves measured at the higher speeds 87 percent and 97 percent 
form a different group. For this group, the phase angles are near 45 degrees 
over the frequency range from 120 to 240 Hz. 

Comparing the transfer function magnitudes calculated with and without 
flow with the measured transfer functions at 33 percent and 40 percent engine 
speed shows good agreement with both theories. However, the transfer function 
calculated for the variable area nozzle with flow fits the phase angle best. 
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Surprisingly, the transfer function calculated without flow is in fair 
agreement with the measurements at all engine speeds. This shows that while 
the nozzle exit impedance model used is adequate at the low flow conditions, 
at the higher flow conditions a more adequate model for the exit impedance and 
more information about the convected entropy is needed, 

CONCLUDING REMARKS 

An acoustic transmission matrix analysis of sound propagation In a 
variable area duct with flow and convected entropy fluctuations was presented. 
The analysis was used with an acoustic impedance and an entropy-pressure 
transfer function exit boundary condition to calculate the acoustic-pressure 
transfer function between the Inlet and the exit of a JT15D core turbofan 
engine nozzle. The one-dimensional model for sound propagation in a duct 
yields good results if the duct exit impedance model is adequate. In spite of 
the scatter in the data, it is possible to see that the impedance model used 
is adequate for low Mach number flows (M less than 0.2). and at frequencies 
from 120 to 720 Hz. It is believed that a more adequate duct exit impedance 
is required at higher Mach numbers and at frequencies less than 120 Hz. 

In order to find the effect of the convected entropy, its magnitude and 
the transfer function between pres« .re and entropy must be known. However, 
for this set of data the effect of convected entropy seems so negligible that 
it is lost in the scatter of the phase angle data at the low engine speeds 
where the exit impedance model used works best. At the higher engine speeds, 
its existence is obscured by the failure of the exit impedance model. These 
results show that; (1) the one-dimensional sound propagation model does apply 
to the 0T15D turbofan nozzle data; (2) in order to study the convected entropy 
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effect In a real engine it will be necessary to measure its magnitude and the 
entropy-pressure transfer function; and (3) the exit impedance model must be 
improved. 
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TABLE I. - TEST CONDITIONS FOR THE SELECTED 0T15D 
NOISE MEASUREMENT POINTS 


Nominal engine 
speed, 
percent of 
maximum 

Core 

Pres- 

C t 1 

Temperature, 

K 

Exit 

speed 

of 

sound, 

m/s 

Air speed, 
m/s 

Mach 

number 

mass 

flow, 

kg/s 

bUi W ) 

atm 

exit 

inlet 

exit 

inlet 

33 

2.3 

0.997 

675 

515 

81 

40 

0.16 

0.08 

40 

2.7 

1.008 

687 

519 

96 

47 

.19 

.09 

55 

3.8 

1.037 

717 

529 

136 

67 

.26 

.13 

87 

6. 6 

1.097 

772 

543 

198 

94 

.37 

.17 

97 

- 

7.8 

1.264 

878 

569 

320 

142 

.56 

0.24 


p 
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Figure 1. - JTISO turbolin engine on test sUnd. 
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Figure 2. • Core acoustic waveguide probe iocations. 
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AXIAL POSITION, cm 
Figure 3. - Core nozzle geometry. 
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Figure 4. - JT15D core nozzle trans- Figure 4, - CcntiRaed. 

mission matrix, T where Y»x-39.5cm> 

“T Y(x-O.C). 
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I(x=U =50.0J/kg K 
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FREQUENCY Hz engine speed and calculated aulo-spectra. 

(a) Magnitude. 

(b) Phase angle. 
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„ (d)JT15D nozzle transfer function 87 

(c) JT15D nozzle transfer function 555^. 

Figure 9. - Continued. 
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